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We present a detailed study of a simple scalar field model that yields non-singular cosmological 
solutions. We study both the qualitative dynamics of the homogeneous and isotropic background and 
the evolution of inhomogeneous linear perturbations. We calculate the spectrum of perturbations 
generated on super-Hubble scales during the collapse phase from initial vacuum fluctuations on small 
scales and then evolve these numerically through the bounce. We show there is a gauge in which 
perturbations remain well-defined and small throughout the bounce, even though perturbations in 
other commonly used gauges become large or ill-defined. We show that the comoving curvature 
perturbation calculated during the collapse phase provides a good estimate of the resulting large 
scale adiabatic perturbation in the expanding phase while the Bardeen metric potential is dominated 
■ by what becomes a decaying mode after the bounce. We show that a power-law collapse phase 

with scale factor proportional to (— t) 2 ^ 3 can yield a scale-invariant spectrum of adiabatic scalar 
perturbations in the expanding phase, but the amplitude of tensor perturbations places important 
CNJ . constraints on the allowed initial conditions. 
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Current cosmological observations, especially measurements of the spectrum of temperature anisotropics and po- 
larisation of the cosmic microwave background, provide strong evidence for the existence of primordial density per- 
turbations on large scales, far larger than the Hubble scale at the time of last-scattering. Observations are consistent 
, with an almost scale- invariant, Gaussian distribution of adiabatic primordial perturbations 

The standard explanation for the origin of the primordial perturbations is that the Universe underwent an accel- 
erated, inflationary, expansion in the early universe |2|. Zero-point vacuum fluctuations in a massless scalar field 
naturally give rise to a scale-invariant distribution of effectively classical field fluctuations on super-Hubble scales dur- 
ing an exponential (de Sitter) expansion. If inflation is driven by a light, slowly-rolling, inflaton field then these field 
fluctuations can be directly related to a comoving curvature perturbation, 1Z, which remains constant for adiabatic 
Oh' perturbations on super-Hubble scales until a given wavelength re-enters the horizon in the subsequent radiation- or 
q , matter-dominated eras. 

It is important to consider whether or not inflation is the only consistent model for the origin of large-scale structure 
in our Universe. An alternative proposal is that the hot big bang and the primordial perturbations originate not from 
an accelerated expansion, but from a preceding collapse. In classical general relativity a collapsing universe may 
be doomed to a singularity, but in a finite quantum theory of gravity, such as string theory, the spectrum of states 
going into the "singularity" would be expected to determine some out-going state. In particular the collapse phase 
could set the initial conditions for a subsequent post-big bang phase. Gaspcrini and Vcneziano originally proposed 
a pre big bang scenario based on the superstring low energy effective action More recently Khoury et al 
proposed the ekpyrotic scenario @ motivated by colliding brane solutions in heterotic M-theory (see also @ for 
later developments of this scenario). Both models invoke a collapsing universe in the 4D Einstein frame where the 
energy density is dominated by minimally coupled scalar fields. 

In a collapsing universe zero-point vacuum fluctuations with a given wavelength arc squeezed outside the rapidly 
shrinking Hubble scale, also naturally giving rise to a spectrum of field fluctuations on super- Hubble scales. However in 
either pre big bang or ekpyrotic scenarios, vacuum fluctuations give rise to a steeply blue-tilted spectrum of comoving 
curvature perturbations during the collapse 0, |t| . If the tilt of the comoving curvature spectrum remains unchanged 
during the transition from pre to post big bang phase then this cannot seed the observed large scale structure in 
our Universe. One would then have to consider alternative sources for structure after a collapse phase, such as the 
curvaton mechanism |Ioj . or consider a different collapse model. Collapse driven by matter or a scalar field with 
pressureless equation of state can produce a scale- invariant spectrum of comoving curvature perturbations 0, 0] . 

Although very general physical arguments have been advanced that the comoving curvature perturbation, 1Z, 
should be conserved on super-horizon scales for adiabatic density perturbations [l3L [HI ITU Uti Il7| there are several 
possible ambiguities that may arise in the transition from collapse to bounce, not present in conventional models of 
an inflationary expansion. At any bounce the comoving Hubble scale diverges and all scales are inside the Hubble 
scale, at least for an instant 0. And the Bardeen metric potential, does not have the same spectral dependence 
as the comoving curvature perturbation during the collapse [13, |2(J, 123, UM > even though the two tilts must 
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coincide for adiabatic perturbations in the subsequent expanding universe. As a result it is impossible for both 1Z and 
$ to obey non-singular evolution equations through the bounce |25| . 

In fact the Bardeen potential on super-Hubble scales grows rapidly in a collapsing phase. This was noted in Ref. @ 
who realised that the divergence of perturbations in the Bardeen's longitudinal (or conformal Newtonian) gauge could 
be eliminated by a gauge-transformation to the spatially flat (or off-diagonal) gauge. However we will see that at 
the moment of bounce the spatially flat gauge itself becomes ill-defined. The comoving curvature perturbation may 
remain small during collapse but it too becomes ill-defined in the vicinity of the bounce in our model. The breakdown 
of perturbation theory in any given gauge does not necessarily signal the breakdown of perturbation theory. It is well 
known that the comoving curvature perturbation becomes singular at stationary values of the field even in single-field 
inflation models, but there are other gauge choices where perturbations remain small |2(j. 

Several authors have studied the evolution of perturbations through a bounce using some analytic form for the 
background scale factor [2?], HE |2t| and/or the evolution equations for perturbations 25]. But this has lead to 
apparently contradictory results depending on which functions are assumed to have a smooth regular form 25] . 

In order to investigate some of these issues associated with the evolution of cosmological perturbations through the 
bounce we will consider a simple scalar field model for a non-singular bounce. We are following the general approach 
envisaged in the pre big bang scenario [3j where higher-order corrections to the effective action will become significant 
at high curvatures or strong coupling and can yield non-singular bounce solutions |3Ct l3lj . In the ekpyrotic or 
cyclic models the transition is supposed to occur at a real singularity in the effective-4D theory. Not only background 
quantities but also perturbations diverge here. So our 4D analysis will not address arguments about the regularisation 
schemes at the singularity proposed so far j^, HE Hi| . 

There have already been a few numerical studies of the evolution of perturbations through non-singular bounce 
solutions including higher-order corrections to the string effective action HE HE HI ■ It can be shown that that 
the null energy condition must be violated for a bounce to occur in general relativity in a spatially flat 4D FRW 
cosmology 1381 . (This is quite different from non-singular classical bounces which have been studied in closed FRW 
models |lllH EE El El El El E3 which do not violate the null energy condition, but require fine-tuned initial 
conditions to obtain a sufficiently long collapse phase 0, H3-) O ne problem with results based on higher-order 
corrections near the bounce is that the evolution close to the bounce becomes dominated by the "corrections" and 
one cannot be sure that even higher-order corrections would not also become important |3lj . Another complication 
is that increasingly higher-order corrections are equivalent to introducing additional fields with non-canonical kinetic 
terms |46L |47| so that the true dynamical degrees of freedom become obscured. 

In order to make the physical degrees of freedom manifest we will work with scalar fields obeying canonical wave 
equations in four-dimensional Einstein gravity. The price we must pay for such simplicity, while still obtaining a 
non-singular bounce, is that one of the scalar fields will be a massless "ghost" field with negative kinetic energy, in 
order to violate the null energy condition near the bounce. Such a ghost field would lead to the presence of serious 
instabilities in any low energy theory |48(. However it does not play any significant role in our model at early or 
late times. Instead it leads to precisely the sort of instability at high energies that we require to achieve a bounce in 
the high-energy regime, destabilising the pre big bang collapse phase and triggering the transition to an expanding 
cosmology. Our scalar field model is intended only as a simple model of a non-singular bounce and not as a realistic 
model of the low energy universe. The early- and late-time (low energy) evolution in our model is dominated by a 
conventional minimally-coupled scalar field with exponential potential. By choosing the dimensionless slope of the 
potential we are able to investigate different power-law solutions for the pre big bang collapse leading to different 
power spectra for the field and metric perturbations A very similar bounce model was studied by Peter and 

Pinto-Neto [22j, but we find very different results. 

In the original pre big bang scenario higher-order corrections to the effective action will become significant at high 
curvatures or strong coupling and can yield non-singular bounce solutions [30l l3l| , at least for some initial conditions 
and some higher-order corrections. In particular it can be shown that the null energy condition must be violated for 
a bounce to occur in a 4D FRW cosmology |38l ]. 

We first study the background evolution of our model using a phase-plane analysis in Section [n] With a particular 
choice of dimensionless phase variables the early or late time solutions correspond to critical points in the phase plane. 
However an alternative choice of variables is better suited to study evolution through the bounce. We identify a class 
of non-singular models which start in a power-law collapse in the asymptotic past, pass smoothly through a bounce 
and approach a power-law expansion. 

We then define linear perturbations about this background model in Section UTTl We present the evolution and 
constraint equations for the perturbations in arbitrary gauge and then identify particular gauge choices that are well- 
suited to study the perturbations at early and late times, or through the bounce, analogously to the different choice 
of background variables used to study the background phase-plane in different regimes. The choice of gauge to take 
the perturbations through the bounce is vital as it must be ensured that the gauge itself is well defined throughout 
the bounce, and we identify problems with gauges commonly used in inflationary models. 
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In section Hvl we evolve first-order evolution equations for the metric and field perturbations and use the constraint 
equations as a consistency check of our calculation. We assume perturbations begin in the quantum vacuum state on 
small scales and use this to calculate the spectrum of perturbations generated on super-Hubble scales during collapse. 
We evolve this numerically through the bounce in order to calculate the resulting spectrum of curvature perturbations. 
We present our results in section in terms of the gauge- invariant variables TZ and <E>. We find that it is possible 
to produce a scale-invariant spectrum of adiabatic perturbations on super- Hubble scales from a collapse phase in our 
simple model. Finally we present our conclusions in Section IV*fl 



II. BACKGROUND 



We consider a spatially flat FRW cosmology containing two minimally scalar fields; one with positive kinetic energy 
and self-interaction potential and one with negative kinetic energy and zero potential energy. The Lagrangian of the 
system is 

£ = -\^^-V{i P ) + \x ll x,»- (i) 

With our choice of signature the field ip has a positive definite energy density and x is a "ghost" field with negative 
energy density. 

We take the potential of the field to be a simple exponential 

V = V cxp(-Akv?) . (2) 

Varying the action with respect to each of the fields and the metric then gives a Klein-Gordon evolution equation for 
each of the fields and an acceleration equation 

/ + 2V + « 2 ^ = 0, (3) 
2h X ' = 0, (4) 

h' = y(V 2 +x' 2 + aV). (5) 



subject to the Friedmann constraint 



where a prime denotes differentiation with respect to conformal time, T). 
We can integrate Eq. (0J to obtain the first integral 



h = T ~ T>X + a V , (6) 



X'=*. (7) 



Thus the negative energy density of the ghost field decreases as the universe expands and grows as it contracts, 
inversely proportional to the square of the volume, x n j a 1 x a 
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The overall equation of state of the two scalar field system is 



P 

w = — 
P 



which can be written as 

3*V 2 - X ' 2 ) 

w = - l + 37^ • (9) 

Thus we have a "phantom" equation of state (w < —1) |43| whenever the total kinetic energy of the fields goes 
negative, i.e., for ip' 2 < x' 2 - 
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A. Asymptotic phase-space variables 



To describe the asymptotic behaviour of the system we define the dimensionless phase-space variables |5C 



Vet. 



nayV 

' s W (1 °> 



VEh ' 

The Friedmann equation JfjJ then gives a constraint equation for the system 

x 2 +y 2 -z 2 = l. (11) 

This reduces the system to become two dimensional, with all the solutions lying on the two-dimensional surface of a 
unit hyperboloid in three-dimensional Euclidean phase space. 

The Klein-Gordon equations © and Q for the two fields and the evolution equation JSJ) give the first-order 
evolution equations for the phase variables 

^ = -Z x {l-x* + z*) + \sJ\y\ (12) 



(/y y(3x 2 -3z 2 -Xd-x), (13) 



dz 

In 



-3z(l - x 2 + z 2 ) . (14) 



where N = In a is the number of e-folds of expansion of the universe. 

This system of equations can be solved numerically to find the phase-trajectories. In figure ^ we show phase 
trajectories for A = 2 projected onto the y-x and z-x planes. Without loss of generality we take x' > 0, as from Eq. Q 
X must either be monotonically increasing or decreasing. 

The fixed points of the system (critical points where ^ = ^ = ^ = 0) are at 



A± : x — x, y — 0, z — ±\/ x 2 — 1 , (15) 
B±: x=j., y = ±y/l-f, z = 0. (16) 

following the notation of [5j . 

^4± are hyperbolic lines of critical points (|x| > 1) where the potential energy V(ip) is negligible. The universe is 
dominated by the kinetic energy of the two fields with ip' oc x' <x h. The scale factor evolves as a power-law with 
respect to conformal time, a cx Ir]] 1 ^ 2 (or a oc It] 1 / 3 in terms of proper time, dt = adrf). 

Critical points B± describe a power-law solution where the density of the phantom field x is negligible, and the 
potential and kinetic energy of the ip field are proportional. This is the well-known scaling solution for exponential 

potentials[Hi! Hi 

acx|i| p , p = 2/\ 2 . (17) 

Note that for a positive potential Vb > (which we assume here) B± only exist for A < ^/6. Thus, the phase-plane 
evolution will fall into two distinct parameter regimes defined by A ^ V6- Within either of these regimes however, 
the precise value of A will not affect the qualitative description of the evolution. Particular cases of interest are A = 2 
for which the universe described by point B± evolves like a radiation-dominated universe with P — p/3, or A = y3 
for which the universe evolves like a dust-dominated universe with P = 0. 

We find that whenever B± exits they are a late-time attractor as JV — > oo, i.e., at late times in an expanding 
universe or early times in a contracting universe. 

We will be interested in the non-zero set of collapsing universes which start at point _B_ in the asymptotic past. 
This solution is unstable as we approach a — > and generic solutions either collapse to a kinetic-dominated singularity 
at A± or bounce, at which point \y\ — » oo. There is another non-zero set of collapsing solutions which start at A + 
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FIG. 1: The phase trajectories in the (x,y) and (x,z) phase-planes for A = 2. 



before bouncing, but we shall focus on models where the energy of the ghost field is negligible, z — > 0, at both early 
and late times. 

The phase-plane H12H14f> gives a useful representation of the evolution at finite h and in particular clearly identify 
the asymptotic fixed points at early and late times. However the variables IjlOfl are liable to diverge at any bounce 
point where h = 0. Thus to study the evolution through the bounce we seek an alternative choice of variables such 
the variables remain finite at the bounce. 



B. Phase-space variables near bounce 



Close to the bounce we will describe the expansion h, and kinetic energy and potential of the if field relative to the 
kinetic energy of the x field. From Eq. J7J we can see that the kinetic energy of the \ field remains non-zero and is 
greatest at the bounce. We can rearrange the Friedmann equation to get 



2a 2 V 6h 2 



/2 

,2 1 ,2 2 /2 ' 

X X « X 



1 = r 72+— —> (18) 

and define phase variables 



<£_ 
X' 



(19) 



P=^, (20) 
X 

7 = ^- (21) 

KX 

which are all non-singular through the bounce. Then Eq. <|18|) gives the constraint 

a 2 +/3 2 - 7 2 = l, (22) 

so all the phase trajectories lie on the two-dimensional surface of a unit hyperboloid. The bounce (7 = 0) corresponds 
to the neck of the hyperboloid (a 2 + f3 2 = 1). 

At the same time it is useful to change the "time" variable to a new variable q defined by 

<m 
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FIG. 2: Phase-space trajectories in the (a, 7) plane for A = 2. Everywhere within the red dashed lines the universe will have 
a phantom equation of state. 



so q is a dimensionless representation of the scalar field x- The advantage of this is that, unlike N, q will be either 
monotonically increasing or decreasing through the bounce and for all time, as can be seen from Eq. Q. 
The system of equations (j2HSJ) then becomes 



da 2 



dq 2 



P , (24) 



T q = ^ ~ — a) ' (25) 

dry 
dq 



3/T . (26) 



Equations and (|26|l can be integrated to give a first integral 

a-^l = C, (27) 

where C is a constant. Thus the solutions evolve along straight lines in the (a, 7) plane, as shown in Figure J5J. 
The only fixed points of this system are at j3 = 0, i.e., they lie on the hyperbolae 

A±: a = ± v / l+7 2 , (28) 

which describes the kinetic-dominated scaling solution, tp' oc x' °^ h. 

As before, the phase plane evolution falls into two categories depending on whether A 5s In terms of the 

variables a and 7 we see that this corresponds to the slope of the trajectories being greater or less that unity. This 
determines whether or not trajectories hit the hyperbolae A±. 

For A 2 < 6 the phase plane looks like that shown in Figure [21 The trajectories fall into different classes according 
to the choice of initial conditions. This is shown schematically in Figure 03 

The evolution of the universe in each of the classes is 

A : power-law contraction (p = 2/A 2 ) — > bounce — > power-law expansion (p = 2/A 2 ) 

B : kinetic dominated contraction — > bounce — > power-law expansion (p = 2/A 2 ) 

C : power-law contraction (p = 2/A 2 ) — > bounce — > kinetic dominated expansion 

D : Big Bang — - > kinetic dominated expansion — ► power-law expansion (p = 2/A 2 ) 

E : power-law contraction (p = 2/A 2 ) — ► kinetic dominated contraction — * Big Crunch 

For A 2 > 6 the phase trajectories evolve as shown in Figure 0| The different classes of solution are shown schemat- 
ically in Figure The trajectories fall into only three classes in this case. The evolution of the universe in each of 
the classes is: 
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FIG. 3: Classification of phase-space trajectories in the (a, 7) plane for A = 2. 




FIG. 4: Phase-space trajectories in the (a, 7) plane for A = 3. Everywhere within the red dashed lines the universe will have 
a phantom equation of state. 



F : kinetic dominated contraction — > bounce — ► kinetic dominated expansion 
G : Big Bang — > kinetic dominated expansion 
H : kinetic dominated contraction — > Big Crunch 
The overall equation of state © of the system is 

2(a 2 - 1) 

w = -l+ V 5 ' ■ 29 

T 

In order to have phantom behaviour (w < —1) we therefore require a 2 < l.For A = 0, the trajectories are a = const, 
and so are straight vertical lines in the 7 — a phase plane. For A ^ no trajectory remains within the phantom region 
and we always have w > — 1 at early times and late times. 

In the second part of the paper we will study inhomogeneous perturbations about the solutions in class A, i.e., 
non-singular evolution for A 2 < 6 where the ghost field is negligible at early and late times. From Eqs. (1221 and 1271) 
we find that this corresponds to solutions with integration constant 



(30) 
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FIG. 5: Classification of phase-space trajectories in the (a, 7) plane for A = 3 



III. LINEAR PERTURBATIONS 



The most general metric of a perturbed FRW universe can be written to first order 



g^u = a 2 in) 



-(1 + 20) 

B b - Sj (1 - 2-0)7^ + 2E Hj + Fi\j + F jH + h t 



B\i — Si 



The scalar quantity <p describes the perturbation in the lapse function, tjj describes the perturbation of the intrinsic 
3-curvature and the scalar shear is given by 



a = E 1 — B . 



(31) 



We need to consider only the scalar and tensor parts (hij) of the metric perturbations, which remain decoupled to 
first-order. Vector perturbations (Si and Fi) can be eliminated by a gauge transformation when the energy-momentum 
tensor only comes from scalar fields [5^| . 

We can also express the perturbed fields in terms of their homogeneous background value plus an inhomogeneous 
perturbation 



<p = Lp + 5ip, 
X = Xo + Sx- 

Under a first-order time-shift, or temporal gauge transformation, 

v = v + e° , 

the scalar metric perturbations transform as [55| 

$ = 4'+ he, 

a = a + C , 
and the first-order field perturbations transform as 

sip = sip-v'e, 

Sx = 5 X - X 'e. 



(32) 
(33) 



(34) 



(35) 
(36) 
(37) 



(38) 
(39) 



The tensor metric perturbations, hij, are gauge-invariant. 
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The linearly perturbed Klein-Gordon equations for the two fields give the wave equations for the field perturbations: 

Sip" + 2h6ip' -V 2 5ip + a 2 V vv 6(p + 2a 2 V v (j)-(p' <j)' -tp' {3ip' ~ V 2 cr) = 0, (40) 
S X " + 2h6 X ' -V 2 5 X -Xo<P' -XoW - V 2 <x) = 0. (41) 

Similarly the Einstein equations can be split into background and perturbed parts, yielding the background Eqs. (J5J 
and jnj), and to first-order 

2 

-3W + </0+V 2 [^ + H = -y(^'o 2 -Xo 2 )-^' + x' 5x' -a 2 V v 6ip), (42) 
W + V = ~(Xo*X-^). ( 43 ) 

2 

(2ti + h 2 )cp + h4>' + r + 2htf = y [ tp'oStp' - x' Sx' ~ Hv'o 2 - Xo 2 ) - a?V, v 5<p ] , (44) 
a' + 2ha + ip - <j> = 0. (45) 

An arbitrary inhomogeneous perturbation can then be decomposed into a superposition of independent Fourier modes 
with comoving wavenumber k, such that V 2 5<p = —k 2 8ip. 

The gauge-transformation equations H35H39|) show that the evolution equations contain an unphysical gauge mode, 
identified with the choice of constant-time hypersurfaces, as well as the physical degrees of freedom. To remove this 
ambiguity one can work with gauge-invariant variables, but there is no unique choice of gauge-invariant variables, 
such as, for example, the curvature perturbation. 

To determine the initial vacuum fluctuations at early times in a collapsing phase and then follow these perturbations 
through the bounce phase we will solve the perturbations in two specific choices of gau ge, e ach adapted to one phase 
of the evolution. Either gauge can be used to define different gauge-invariant variables |5g . Because we have gauge- 
invariant definitions in each phase we can consistently follow the perturbations through from early to late times. 

Ultimately we will present our final results in terms of the two most commonly used gauge-invariant curvature 
perturbations: 

n = ^ h ^' 5 ri 5x \ (46) 

iff* - x l 

$ = ^j + ha. (47) 

Equation (|4f>|) is the natural generalisation of the usual comoving curvature perturbation for two scalar fields |57ll58| to 
the case where one of them is a ghost field, while Eq. H47|) is the curvature perturbation in the longitudinal gauge |55| . 
also called the Bardeen potential. 



A. Uniform curvature gauge 



In order to describe the early time behaviour close to the fixed point we will first work in the uniform curvature 
(or spatially flat) gauge in which the metric perturbation, ij) vanishes. This is the gauge commonly used to describe 
vacuum fluctuations generated during inflation |59l l60l l61| . 

From Eq. we see that -0 = corresponds to a specific time-shift from an arbitrary gauge 

e = (48) 

This enables us to give a gauge-invariant definitions of the remaining metric perturbations |56| and, in particular, the 
scalar field perturbations in the uniform-curvature gauge 

8tp\* = 6<p + fy, (49) 
<f X |v, = Sx + jiP. (50) 



In the rest of this sub-section we will drop the tilde and "0-subscript, but it should be remembered that all perturbations 
correspond to those in the uniform-curvature gauge. 



10 



Using Eas. (|42H45|l to eliminate <fi and tp, the perturbed Klein-Gordon equations in the uniform-curvature gauge can 
be written as two coupled second-order equations for the field perturbations: 



Sip" - V 2 Sip + 2hSip' + 



a 2 V w + k (p f Q ( - 2) - 2^ 



Sip 



•>/ v, ' -,■'-/•> i ^' \x-. k2 Xo h 1 ' i\ x 



Sx" - V 2 <5 X + 2h5 X ' - k 2 x ' 2 (2 + -^)S X = 



(51) 
(52) 



These evolution equations in the uniform-curvature gauge have a simple form which is convenient for studying the 
evolution during monotonic expansion or collapse. In particular we will use this gauge to describe the initial quantum 
fluctuations during a collapse phase. 

However these equations are not suitable for evolving the perturbations through a bounce. The evolution equa- 
tions l|51ll and H52fl contain singular coefficients as h — ► 0. This is not surprising as from Eq. H48|) we see that the 
uniform-curvature gauge is itself ill-defined at a bounce as h — ► 0. Instead, we have to find another gauge which 
remains well-defined at the bounce in order to evolve the initial vacuum perturbations through to an expanding 
phase. 



B. Uniform x-field gauge 

Analogously to the background solutions where we used the x _ field as a useful time coordinate to study evolution 
through the bounce, we will use the uniform- x gauge to follow linear perturbations through the bounce. 
From Eq. (|39|l we see that to set S\ = we must choose 

e = H. (53) 

X 

The remaining t^-field perturbation in the uniform- x gauge has the gauge invariant definition 

S~ip\ x = Sip-^S X . (54) 
X 

The scalar metric perturbations cj>, ip and a in the uniform- x gauge also have gauge-invariant definitions. In particular 
the curvature perturbation on uniform- x hypersurfaces is given by 

$\x = $ + n s x- ( 55 ) 

X 

In the rest of this sub-section we will drop the tilde and x-subscript, but it should be remembered that all perturbations 
correspond to those in the uniform-x gauge. 

The perturbed Klein-Gordon and Einstein equations l|40H45|) . then give us the four evolution equations for the four 
perturbation variables dip, tp, </> and a in the uniform x-field gauge 

6fc 2 

K5<p qq + 2 k ^ + 3/3 2 A 2 k^-6A/3 2 (/) = 0, (56) 

k 2 x' 

ip m - ^LnSip, q - 3 ^f nSip + j4> tq + 3f3 2 4> = 0, (57) 
vo 2 

6fc 2 

<f>,«+3t« + ^2 s = 0. (58) 
« X 

s tq +4js + ip - 4> = 0, (59) 

plus two constraint equations 

l<t> + 1>, q -^K5ip = 0, (60) 
(a 2 - 1)0 - ^=KS<p, q + - - W> + 7«) = . (61) 
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Here we have (as for the background evolution) we have written time derivatives with respect to the dimensionless, 
and monotonic, variable q defined in Eq. 1|23[) . We have written background quantities in terms of the dimensionless 
variables a, (3 and 7 defined in Eqs. 

We have rescaled the shear perturbation variable a in terms of 

s = -a, (62) 

7 

so that all gradient terms enter through terms of order (k/ kx') 2 which is the square of 7 times the ratio of the Hubble 
length to the wavelength. Thus we can obtain solutions order-by-order in terms of (k/nx') 2 to define a long-wavelength 
limit. Although the Hubble length diverges at the bounce we have from Eq. that (k/Kx') 2 oc a 4 . This remains 
finite and actually reaches a minimum at the bounce. Thus we find that the long-wavelength solution, neglecting 
terms of order (fc/ kx') 2 is an excellent solution close to the bounce for scales of interest. 

In principle we could use the constraint equations ifBTljl and l|6*T|) to eliminate cj) and s respectively and leave two 
second-order equations for Sip and tp, similar to what was done for the uniform curvature gauge perturbations in the 
preceding subsection. However as 7 — > approaching a bounce this approach would be very sensitive to any small 
numerical error. We will not eliminate cj> and s and instead will use their first-order evolution equations to evolve them 
numerically through the bounce. The constraint equations (|6Ull and l|61l) then provide a useful consistency check for 
our numerical integrations. A similar argument in favour of using coupled first-order equations to study perturbations 
through a bounce was recently given by Cartier [3(j . 



C. Tensor perturbations 



The tensor part of the perturbed metric, hij, is transverse and trace-free and is invariant under gauge transforma- 
tions. We expand hij in plane-waves 

hij (x, V ) = (2n)- 3 f d 3 k Y, 5 9l{ri) % (k, x) , (63) 

p=X,+ 

where e?. is the polarisation tensor and p = x, + represents the two independent polarisation states. Sg£' + is the 
scalar amplitude of each state for wave-mode k. For both polarisation states this scalar amplitude obeys the simple 
wave equation for a free scalar field in an FRW spacetime, 

Sg" + 2h8g' + k 2 Sg = , (64) 

where k 2 = k.k. Here and from now on the subscript k and the superscript x , + are assumed for Sg. 

Under transformation to derivatives with respect to our monotonic time variable q defined in Eq. (|23|l . the tensor 
evolution Eq. I|64|) becomes 

6fc 2 

s 9,qq + a S 9 = • ( 65 ) 

K Z X 

This is the form of the equation that we use to evolve the tensor perturbations numerically through the bounce. 



IV. VACUUM FLUCTUATIONS 



A. Early-time behaviour in uniform-curvature gauge 



For any power-law solution described by a critical point in the dimensionless phase-space we have (p' oc h and the 
right-hand-sides of Eqs. I|51() and l|52() vanish. Thus the field perturbations Sip and 6x m the uniform-curvature gauge 
decouple in the early-time limit where the background solution is described by Eq. 1)170. In this case the evolution 
equations for the field perturbations can be written very simply in terms of the canonical Sasaki/Mukhanov variables 



= aSip , 
= aSx, 



(66) 
(67) 
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For each Fourier mode with comoving wavenumber k, Eqs. i|51|) and (|52|l can then be written as 

k 2 - — ]ui = 0, (68) 



where i = p,x- The evolution equation for the tensor perturbations Eq. 16411 can be written in exactly the same form, 
where 



aSg 
2k 



*9 ~ — ■ ( 69 ) 



For power law expansion or collapse described by Eq. I|17|l the equation 1|68[) has the general solution |5^| 

u = y^NK^ON) + u-H^iM)] , (70) 



where the index 

3 1 



p-1 



(71) 



Although the scalar fields ip and x are rea l scalar fields, we have chosen to Fourier expand the perturbations in terms 
of complex Fourier modes oc e ±lkx . Thus we work with the complex modes functions <|TU|) which must obey the 
condition u{— k) = u*(k). 

At early times (r) <C — fc _1 ) we have free oscillations u oc e ±tkri . The simplest assumption is that the pertur- 
bation variables, u% t 2, start in the flat space-time vacuum state with only positive frequency modes, satisfying the 
normalisation condition 



Thus at early times (krj — ► — oo) we require 



u* u = i . (72) 



! e- jfc " . (73) 



(2k) 

Using the vacuum solution (|73(l at early times to set the initial conditions for the general solution (|70|l yields 



\k^H\l (\k V \). (74) 



We use this analytic vacuum solution for the three perturbation variables to set up the initial conditions for the 
perturbations in our numerical code. We start the code far away from the bounce when we are still well within the 
power-law regime (p ~ 2/X 2 ). However, as we saw in Section^ this power-law solution is unstable and as the universe 
collapses the evolution will eventually diverge away from the critical point. 



1. Initial power spectra 

Sip, Sx and 5g are independent Gaussian random fields whose distribution will be completely determined by the 
power spectra, defined for a scalar Sx as ; 68| 

V. 

The tilt of the spectrum is defined by 

An 

so that An x — for a scale-invariant spectrum. Note that the conventional spectral index for scalar metric perturba- 
tions is given by n = 1 + An x . 

As Sip and Sx are independent fields in the initial power-law regime, we can find the power spectrum for any 
quantity by evolving separately different modes corresponding to an initial perturbation in ip or x- The power due 



k 3 \Sa 



2vr 2 



= dln(V x ) 
~ dln(fc) 



(75) 



(76) 
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to each mode can be summed to obtain the final spectrum, tp-modes correspond to initial vacuum perturbations of 
the (/3-field on the uniform curvature hypersurface with the x perturbation (and its derivative) initially zero on this 
hypersurface. x-modes correspond to the opposite situation, where there is an initial perturbation of the x-field but 
no perturbation of the tp field (or its derivative) on the uniform curvature hypersurface. 

Any complex variable can be split into two parts whose phase differs by 90°, which typically might be real and 
imaginary parts. We find it convenient to split the Hankel function, hPa, found in our general solution for u, Eq. 174fl . 
into real Bessel functions of the first and second kind: 

H^(\kri\)=J {vl (\h 1 \)+iY {vl (\h 1 \). (77) 

We will thus refer to these as the J and Y modes. From the definition of the power spectrum Eq. |Z5t it can then be 
seen that 

V x = P XJ + V XY (78) 

Thus, the J and Y parts can be evolved separately to determine the total power spectrum. This gives us four 
independent 'modes' or sets of initial conditions (<pj, fYi Xj? an d Xy) f° r which we present our numerical results. 

The Bessel functions of first and second kind describe regular and singular solutions on large scales (\kr)\ — > 0) 
behaving as 



I (\krj\_\ 



M 



V 2 

Thus, from Eqs 174|) . we see that the power spectrum from vacuum fluctuations during collapse on large scales yields 

V u = C3(M)fc 2 |fery| 1+2 M + C Y {\v\)k 2 \kr,\ l - 2 M , (81) 

where 

CAM) s 2M + l^r(M + D - (82) 

7T2 

The spectral tilts for each term in the spectrum l|81|l are then [Til |59| 

An UJ = 3 + 2|i/|, (84) 
An UY = 3-2|i/|. (85) 

In a contracting universe the second Y-term in Eq. (|81|l should quickly dominate on super-Hubble scales as kr\ — > 0. 
However, we will evolve both modes numerically to see whether this will remains true through the bounce or whether 
the less dominant mode in the collapsing phase will come to dominate in the expanding phase. 

We will present detailed numerical solutions for interesting physical choices of the scalar field slope, A, and hence 
power-law p in the collapse phase that give simple analytic solutions for the perturbations during the initial power-law 
collapse phase. 



2. X — 2: Power-law solution with p = 1/2 

Taking A = 2 means that the scale factor will follow the same evolution as a radiation dominated universe with 
W = 1/3 and p = 1/2 at the critical points B±. 

The Hankel function index in Eq. (|71|l is v = —1/2 and for this special case the first Hankel function has the simple 
form 

j o 

Hr{z) = \ — (sin(z)-icos(z)). (86) 
2 \ irz 
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and the vacuum solutions (|74|l become 



u=-^- T (cos(\kr]\)+ism(\kr]\)) (87) 
(2k) 2 



In this case the l"-part of the solution corresponds to the real part of the solution which will dominate in the collapsing 
universe as kr\ — > 0. 

The spectral tilts of the Y and J parts of u are then given by 



An UY = 2, 

An Uj = 4. (89) 

Both modes have steep blue tilts. 



3. A = \/3-' Power-law solution with p — 2/3 



Taking A = V3 means that the scale factor will follow the same evolution as a matter dominated universe at the 
critical points B±, with w = and p = 2/3. 

The Hankel function index in Eq. (|71|l is v = —3/2 and for this special case the first Hankel function has the simple 
form 



H i%( z ) = ( 1 - 72 ) ( cos ( z ) + isin(z)) . 



and the solution for u will be 



(2fc)3 



cos(|fc?7|) — 



sin(|/c77|) 
1**71 



1(1*7,1) 



COS 



(\ky\) \ 



M J 



(90) 



(91) 



This is the same mode function as is found for a vacuum fluctuations of a massless field in de Sitter inflation |llj |. 
Here the Y"-mode is the imaginary part of the solution will dominate in a collapsing universe as krj — > 0. 
The spectral tilts of the Y and J parts of the field fluctuations on large scale are 

An UY = 0, (92) 
An Uj = 2. (93) 

The J part has a steep blue tilt, but the Y part, which dominates as krj 0, has a scale-invariant spectrum [TH Il2l l62l | . 



B. Numerical evolution through the bounce 

The initial conditions for the background variables are set so that the evolution follows a trajectory in the phase- 
plane that lies within region A, as defined in section (|IIB(l . Thus the evolution starts close to the critical point B_ 
where it follows a power-law collapse solution. It eventually diverges away from this critical point as the universe 
collapses and the (negative) kinetic energy of the \ field grows than that of the <p field. The scale factor smoothly 
evolves through the bounce and finally moves towards the critical point B + in the expanding phase, where once again 
the model approaches a power-law solution with p = 2/ A 2 . 

We use Eqs. (I56I59|) to evolve the scalar perturbations in the uniform x-field gauge through the bounce using q 
as our monotonically increasing time variable. This requires initial conditions for Sip, Sip < q , ip, ip >q , 4> and s in each 
mode. The first four initial conditions are found from our analytical solutions in section l|IV All . However, as these 
solutions are found in the uniform curvature gauge it is necessary to transform them to the uniform x-ficld gauge. It 
is also necessary to change the time variable from r/ to q. The appropriate gauge transformation Ij53(l gives 

Sip\ x = 8ip\^ - aSx\i> , (94) 
S<P,g\ x = -^7r-&xU + —A S( P'U ~ a5 xU) , (95) 
V| x = ^=*XU, (96) 

V, ? | x = ^SxU + ^S X %. (97) 

V D X 
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The other two initial conditions required, for <f> and s, are found from the constraints (|60I61|I . Thereafter the constraint 
equations serve as a consistency check of our numerical solutions. 

The dimensionless metric perturbations ip and (j> remain small in the uniform x-field gauge and well-behaved through 
the bounce. Typical numerical solutions for the curvature perturbation tp is shown in Figure[S]for A = 2 and FigureEl 
for A = V3. The perturbed lapse function (f> closely follows the curvature perturbation on large scales, where we find 
4> + 3ip ~ constant, as expected from Eq. (|58fl. 

The shear potential a does grow rapidly in the collapse phase, reflecting the instability with respect to shear of the 
power-law solution. From Eq. (|44|) we see that in any gauge where 4> and ip remain small, we find a oc a~ 2 as the 
collapse proceeds. However the shear only enters the evolution equation for (f> through its spatial divergence and 
this term remains small as the bounce is approached. We have checked that the magnitude of the shear term in, for 
example, the energy constraint equation (|61[) . is always much smaller than, for instance, the energy density of either 
of the scalar fields. Note that the total density, and hence the isotropic expansion rate, goes to zero at the bounce, 
so the anisotropic shear perturbation cannot be smaller than the isotropic expansion/collapse at this point but it is 
always small relative to the largest terms in the equation so we believe our perturbative analysis is justified. 

All the fc-dependent terms in the evolution equations l|56H59|l in the uniform-^ gauge become small as we approach 
the bounce and we find that the k = solution is an excellent approximation through the bounce for modes which 
leave the Hubble scale well before the bounce (k <§; (aH) max ). Note the difference with models of a classical bounce 
in closed FRW geometries, where finite scale effects may be significant as modes re-enter the horizon during the 
bounce [l8|. 

At late times, in the expanding phase after the bounce, the curvature perturbation in the uniform x-field gauge 
gradually grows and eventually may become larger than one. This is because x' — > as we asymptotically approach 
the power-law solution once again and the gauge shift l|53|) becomes large. Thus at late times in the expanding 
universe we should transform back to some well-behaved gauge, such as the longitudinal or uniform-curvature gauges, 
where the metric perturbations remain small. Nonetheless the uniform ^-field gauge does the job of evolving the 
initial vacuum fluctuations from the collapse phase smoothly through the bounce into an expanding phase where the 
large-scale behaviour is unambiguous. 

The constraint equations (|6U|) and l|61ll are satisfied to one part in ten thousand or better for all our numerical runs 
(see Figures ITT1 and IT5|) . 

V. RESULTS FOR GAUGE-INVARIANT VARIABLES 

In order to compare our results with previous discussions of the expected evolution of perturbations through a 
cosmological bounce, we reconstruct the gauge-invariant variables TZ and <& defined in Eqs. (|4(j[) and (|47|l from the 
variables evolved numerically in the uniform- x gauge. 

For our two chosen values of A = 2 and A = v3 f we find the power spectra of these two gauge-invariant variables 
for each mode when the background is asymptotically evolving as the power law collapse and expansion. We compare 
the spectral tilts of the power spectra before and after the bounce and we see whether the mode dominant before the 
bounce remains so after it. 

Tensor perturbations 8g are automatically gauge-invariant and provide an interesting comparison with the scalar 
metric perturbations. Tensors also turn out to place an important constraint on models producing a scale- invariant 
spectrum of large-scale perturbations. 

A. Comoving curvature perturbation 

Typical evolution of the comoving curvature perturbation, TZ [defined in Eq. I|4fi[) ]. through the bounce is shown in 
Figures[7|for A = 2 andllUlfor A = The four lines correspond to the four modes identified in section llVAl TZ only 
remains constant on large scales for the (p,j mode when A = 2, when this corresponds to an adiabatic perturbation 
of the power-law solution. All the other modes represent non-adiabatic perturbations about the initial power-law 
solution on large scales and have non-constant TZ. In particular the dominant mode, <py, describes non-adiabatic 
perturbations of the ip field. Perturbations in the x field are always sub-dominant in terms of their contribution to 
the comoving curvature perturbation. 

During the bounce phase we see two singular points in the evolution of TZ. This occurs whenever the total kinetic 
energy vanishes, ip 12 — x' 2 = (corresponding to a 2 = 1 in the phase-plane), as can be seen from the definition of 
TZ given in Eq. Ij46(l . At these points the total four- velocity of the two-field system is momentarily ill-defined. This 
makes the comoving gauge unsuitable for evolving the perturbations close to the bounce. We emphasise though that 
we can nonetheless trust the evolution of TZ shown in Figures J7| and (|10|) because here the comoving curvature has 
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been constructed from perturbations in the uniform-^ gauge which remain well-behaved throughout the bounce, i.e., 
we do not require that the perturbations are necessarily small or well-defined in any other gauge. 

After the bounce the comoving curvature perturbation rapidly settles down to a constant value for all modes on 
all super-Hubble scales, representing the resultant adiabatic curvature perturbation on large scales. Thus all the 
modes contribute a finite amplitude to the large-scale curvature perturbation in the expanding phase, and this is 
approximately equal to the amplitude of the comoving curvature perturbation at the end of the collapse phase (but 
not in general the amplitude at Hubble-crossing during collapse). The dominant mode at late times is the ipy mode 
that also dominates the comoving curvature perturbation in the collapse phase. 

Finally the perturbations "re-enter the horizon" in the expanding universe and oscillate for k > aH. 

B. Bardeen potential 

The typical evolution of the Bardeen potential, <& defined in Eq. (|47|l . which describes the curvature perturbation 
in the longitudinal gauge, is shown in Figures Elfor A = 2 andCUfor A = \/3. 

All modes show the same rapid growth of <& during the collapse phase. We can interpret this as being due to the 
rapid growth of the shear in the uniform-^ gauge, a\ x oc a -2 , while the curvature perturbation ip\ x remains small. 
According to the definition of $ in Eq. (|47fl this yields a large curvature perturbation in the longitudinal gauge 

$ oc A ■ (98) 

Our numerical evolution for A = 2 and A = v3 during the power-law collapse phase confirms this, yielding $ oc |t?|~ 3 
and $ oc |?y| -5 respectively. 

Thus the curvature perturbation in the longitudinal gauge becomes large in a collapse phase meaning that this 
gauge too may not be reliable gauge for performing a perturbative calculation through a cosmological bounce. Once 
again though we emphasise that our calculations of the Bardeen potential are reliable as we reconstruct the Bardeen 
potential from perturbative calculations in the uniform-^ gauge which remain small through the bounce. 

After the bounce the dominant modes contributing to the Bardeen potential, ipy and xy, remain dominated by the 
shear, a\ x , which now decays in the expanding phase. Eventually the Bardeen potential settles down to a constant 
value on super-Hubble scales, before finally oscillating after re-entering the Hubble scale. For A = V3 the sub-dominant 
modes, tpj and xj, have much smaller shear (indeed $ stays small for these modes right up until the bounce) and 
these modes rapidly settle down to a constant value soon after the bounce, but this final amplitude is smaller than 
that due to the dominant ipy and xy modes. 

In summary we find that the amplitudes of the Bardeen potential during collapse are not a good predictor of the 
amplitude of the eventual adiabatic curvature perturbation on large scales after the bounce, even though it becomes 
much larger than the comoving curvature perturbation during the collapse. It is the growth of the shear that dominates 
the behaviour of $ during the collapse and this shear then decays in the expanding phase. Only once the shear has 
decayed, sometime after the bounce, does $ settle down to a constant value on super-Hubble scales in the expanding 
era. 

For the case of A = 2 corresponding to a radiation-like equation of state w = 1/3, we find that once the Bardeen 
potential has settled down to a constant value on super-Hubble scales, we have $ = (2/3)7?.. Similarly for A = \/3, 
corresponding to an equation of state w = 0, we find $ = (3/5)7?, as expected for the growing mode of adiabatic 
density perturbations on large scales. 

C. Spectral tilts 

1. \ — 2: Asymptotic power-law solution with p = 1/2 

Table H] shows the initial spectral tilts of the power spectra of 1Z and $, during power-law collapse with A = 2. 
All the modes of 1Z have a blue spectrum, but the dominant mode of $ has a red spectrum. Table [I] also shows 
the spectral tilts of the super-Hubble power spectra of TZ and 4> when they have settled to constant values after the 
bounce. The tilts of the spectra of TZ for each mode remain the same after the bounce as they were initially. The 
spectral tilt of $ for every mode however becomes the same as that for the corresponding 1Z spectra 



(99) 
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FIG. 6: Evolution of ip\ x shown as a function of expansion, N = ±ln(a/ao), for bounce solution with A = 2 model for each 
of the four independent modes: <py (black solid line), tpj (red dotted line), \y (blue dashed line), and \j (purple dot-dashed 
line) . 



mode 


An Kin 


An$ in 


An7j out 


An$ ollt 


ip Y 


2 


-2 


2 


2 


<pj 


4 





4 


4 


XY 


2 


-2 


2 


2 


Xj 


4 





4 


4 



TABLE I: Initial spectral tilts for different modes during collapse and final tilts after bounce for model with A = 2. 



Our model with A = 2 is very similar to the bounce model studied in Ref. 27], but our results are strikingly 
different. The authors in that paper considered a radiation fluid which shows the same initial red spectrum for $. 
But they found that the final spectrum for <!> retained the same red spectrum for $. There are subtle differences 
between a radiation fluid and a A = 2 scalar field. In particular the bounce must be symmetrical with perfect fluids, 
whereas our bounce is only symmetrical for C — and our scalar field supports non-adiabatic pressure perturbations. 
But the main difference may come from the fact that a second-order evolution equation for $ was used to follow the 
evolution of the metric perturbations in Ref. . Perturbations in $ become large and we have instead followed the 
evolution of perturbations in a gauge which remains well-behaved and then reconstructed the value of 4>. We have 
not used the constraint equations to eliminate variables and have instead kept used these as a consistency check of 
our numerical integration. 



2. A = %/3: Asymptotic power-law solution with p = 2/3 

Table ILTI shows the initial spectral tilts of the super- horizon power spectra of TZ and <£>, during power-law collapse. 
The dominant mode of TZ has a scale-invariant spectrum |llL Il2j ] , whereas the spectrum of the dominant mode of <f> 
is red. Table [H] also shows the spectral tilts of TZ and $ when they have settled to constant values on super-Hubble 
scales after the bounce. As for A = 2, the tilts of the spectra of TZ for each mode remain the same after the bounce as 
they were initially. In particular the dominant TZ mode retains its scale-invariant spectrum. The spectra of <E> again 
become the same as TZ after the bounce on large scales 



An$ out = An$ m + 4 = Ann 



(100) 
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FIG. 7: Comoving curvature perturbation 1Z calculated from numerical solution for each of the four modes shown in FigureEl 
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FIG. 8: Bardeen metric potential $ calculated from numerical solution for each of the four modes shown in Figure |S| 



Thus we have shown that a power-law collapse phase with p = 2/3 can yield a scale- invariant spectrum of curvature 
perturbations on super-Hubble scales in the subsequent expanding phase. 



D. Tensor perturbations 

The smooth evolution through the bounce of the amplitude of the tensor metric perturbations is shown in Figure lT2l 
Tensor perturbations grow rapidly during the collapse phase but then rapidly settle down to a constant value after 
the bounce. 
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FIG. 9: Evolution of tp\ x shown as a function of expansion, N = ±ln(a/ao), for bounce solution with A = vo model for each 
of the four independent modes: ipy (black solid line), tpj (red dotted line), \y (blue dashed line), and \J (purple dot-dashed 
line) . 
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FIG. 10: Comoving curvature perturbation 1Z calculated from numerical solution for each of the four modes shown in Figure^] 
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TABLE II: Initial spectral tilts of different modes during collapse and final tilts after bounce for model with A = v3. 
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FIG. 11: Bardeen metric potential $ calculated from numerical solution for each of the four modes shown in Figure 



The spectral index of the tensor perturbations is the same as that of the comoving curvature perturbation 1Z during 
collapse and, like the comoving curvature perturbation, the spectral index is unaffected by the bounce. 

During power-law collapse there is a simple relation between the amplitude of tensor perturbations and the comoving 
curvature perturbation: 



M 

Vn? I 

/ collapse 



8k 2 



For collapse solutions described by point B, we have, from Eqs. i|16|l and Q17[l. 

collapse P 



(101) 



(102) 



In the model with A = ^/i which yields a scale-invariant spectrum of both scalar and tensor perturbations, this 
tensor-scalar ratio is 24 during collapse. 

However both the comoving curvature perturbation and tensor perturbation receive a scale-independent amplifi- 
cation during the bounce. This amplification is dependent upon the details of the bounce, and in particular the 
parameter C defined in Eq. I|27|) which identifies the phase-space trajectory. We show in Figure ^5] the resulting 
tensor-scalar ratio after the bounce. The tensor-scalar ratio is maximised for a symmetric bounce, C = 0, which 
leaves the tensor-scalar ratio (|102|) unchanged. For an asymmetric bounce, the scalars are more strongly amplified 
than the tensor modes and the tesnor-scalar ratio becomes small for the most asymmetric bounces. 

Observations require that the primordial tensor-scalar ratio is less than 0.81 [l|. We see that this places a severe 
constraint on the allowed parameter |C| > 0.67, given that from Eq. I|30[) we require |C| < 0.71 for the collapse 
phase to end in a non-singular bounce. Although we have a tight constraint in our particular bounce model, it is 
not clear how problematic this might be in the more general case as asymmetry through the bounce can decrease the 
tensor-scalar ratio. 



VI. CONCLUSIONS 



We have presented a simple model of a four-dimensional FRW universe that describes non-singular bouncing 
solutions. To do so within the context general relativity it is well-known that the null-energy condition must be 
violated. To achieve this we have introduced a massless "ghost" field with negative kinetic energy. A similar role 
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FIG. 12: Evolution of Sg shown as a function of expansion, N — ±ln(a), for a bounce solution in the A = \/3 model for each 
of the two independent tensor modes: gy (solid black line) and gj (dashed red line). 
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FIG. 13: Tensor-scalar ratio V g /V-R, after the bounce shown against dimensionless parameter C defined in Eq. 1271 . 

The red dashed line shows the observational limit. 



is played by higher-order corrections to the low-energy string effective action which have recently been studied by 
Cartier |36| who found qualitatively similar results. In order to keep the dynamics as transparent as possible we 
have restricted ourselves to a canonical (second-order) scalar field Lagrangian albeit with the wrong sign. Like the 
higher-order loop corrections, our ghost field only plays a significant dynamical role in the high energy regime close 
to the bounce. 

We have performed a phase-plane analysis of the qualitative behaviour of expanding and contracting FRW models 
with both a canonical scalar field (with positive energy density) and an exponential potential plus the massless ghost 
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FIG. 14: The violation of constraint equations 16011 and 1611 shown as a function of expansion, TV = ±ln(a), for a bounce 
solution in the A = 2 model for each of the four independent scalar modes: tpy (black solid line) , tp j (red dotted line) , \y (blue 
dashed line), \J (purple dot-dashed line). The y-axis of the left hand graph shows the logarithm of the absolute value of the 
left hand side of Eq. 16011 divided by the sum of the absolute values of each term in the equation. The y-axis of the right hand 
graph shows the same thing for Eq. \F>1}. 




FIG. 15: The violation of constraint equations l|60[l and ilO 1 1 shown as a function of expansion, N = ±ln(a), for a bounce 
solution in the A = v3 model for each of the four independent scalar modes: ipy (black solid line), tpj (red dotted line), xy 
(blue dashed line), \J (purple dot-dashed line). The y-axis of the left hand graph shows the logarithm of the absolute value of 
the left hand side of Eq. 1)600 divided by the sum of the absolute values of each term in the equation. The y-axis of the right 
hand graph shows the same thing for Eq. Hilt . 



field. Away from the bounce we have usual power-law expansion or contraction for a universe dominated by a scalar 
field with exponential potential |52l l53l l54|. The dimensionless slope of the potential, A, determines the asymptotic 
power-law index p = 2/A 2 and the equation of state w = (A 2 — 3)/3 for sufficiently flat potentials with A 2 < 6. 

It is known that for A 2 < 6 these power-law collapse solutions describe the generic early time behaviour but are 
unstable at late times to kinetic-dominated collapse with w = 1 [Hll l63| or to shear anisotropy [64| . Some kind of 
instability is a necessary to allow an escape from the singular power-law collapse. For steep potentials with A 2 > 6 we 
find no non-singular solutions. But for flat potentials with A 2 < 6 we find a class of non-singular solutions that begin 
in the asymptotic past as a power-law collapse, approach a minimum value of the scale factor and then approach a 
power-law expansion at late times. 

This class of non-singular solutions provides a simple bounce model in which we then study the evolution of linear 
perturbations. We assume that these perturbations begin in the flat-spacetime vacuum state at early times and use 
this to predict the amplitude of large scale (super-Hubble) adiabatic density perturbations after the bounce. 

Although both the Bardeen potential $ and comoving curvature perturbation 1Z evolve smoothly at the instanta- 
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neous bounce, they both have potential problems. The Bardeen potential grows rapidly during the collapse phase and 
is much greater than unity for the modes of interest, calling into question the validity of linear perturbation theory in 
the longitudinal gauge [23, E3 • On the other hand the comoving gauge becomes singular whenever the total kinetic 
energy of the scalar fields vanishes (2(| , which occurs shortly before and after the bounce in our model. 

We can however find a gauge in which the metric perturbations remain small and linear perturbation theory 
should be valid. In our model the massless ghost field is monotonic and provides a suitable time variable, defining 
a well-defined choice of gauge (the uniform field gauge) for perturbations. We note that we solve the regularized 
unconstrained evolution equations for the perturbations, using the constraint equations as a check of our numerical 
integration. This avoids potential problems that could arise from using the constraint equations to eliminate variables 
whose coefficients vanish at the bounce [3(|. We can then reconstruct the usual gauge- invariant variables from 
calculations in this gauge. 

We find that the comoving curvature perturbation calculated during the collapse phase is a good estimator of the 
resulting adiabatic density perturbation on super-Hubble scales in the expanding phase. Note that the value at the end 
of the collapse phase is not in general the same as the that at Hubble-crossing during collapse (unlike the usual case in 
an inflationary expansion). This is because the dominant mode during collapse is a non-adiabatic field perturbation. 
However the comoving curvature perturbation remains constant on super-Hubble scales after the bounce. 

By contrast the Bardeen potential calculated in the collapse phase is not a good indicator of the subsequent large 
scale adiabatic perturbation in the expanding phase. It is dominated by the shear in the uniform x-ficlcl gauge, which 
grows rapidly in the collapse phase and then decays away after the bounce, until the Bardeen potential settles down 
to a constant value on super-Hubble scales, related to the comoving curvature perturbation by the usual expression 
$ = 3(1 + w)lZ/(5 + 3w). This basic point has been made by many authors |9L llTt fl9L |21| and here we demonstrate 
it explicitly for our simple bounce model. 

In particular we find that the spectral tilt of the final adiabatic curvature perturbation on super-Hubble scales in 
the expanding phase is given by the spectral tilt of the initial comoving curvature perturbations (or equivalently the 
scalar field fluctuations in the spatially flat gauge). We therefore have given an explicit model for a power-law collapse 
with A = V3 which yields a scale-invariant of adiabatic perturbations on super-Hubble scales in the expanding phase. 

The amplitude of the comoving curvature perturbation on the smallest scale (that is the minimum value of the 
Hubble scale during collapse) is set by the energy scale relative to the Planck scale (where vacuum fluctuations in the 
metric become of order unity) 

V K (k = aff raax ) ~ (jj^j • (103) 

This is very similar to the amplitude of gravitational waves (tensor perturbations) produced during collapse. The 
tensor-scalar ratio is an important constraint on conventional inflation models and turns out to be a severe constraint 
on our collapse model when A = a/3- 

Pre big bang or collapse models have been criticised for the degree of fine-tuning required to set up the initial 
conditions at the start of any collapse phase [63l Irjri l67| . The dimensionless parameter C, defined in Eq. (|27() , is a 
first integral of the equations of motion which selects the trajectory in the homogeneous phase-space. Although any 
value for |C| < 1/V2 yields a non-singular bounce for models with A = a/3, only trajectories with |C| > 0.67 give 
an acceptably small tensor-scalar ratio. This provides a quantitative measure of the degree of fine-tuning of initial 
conditions required to get an acceptable primordial density perturbation within the class of models considered in this 
paper. 
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